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Abstract 

ft 1 

We present the QCD corrections to the matching conditions of the AB = 1 magnetic 

Qh' 

and chromo-magnetic operators in the Standard Model and in two-Higgs doublet 

(— ! , 

models. We use an off-shell matching procedure which allows us to perform the 
computation using Taylor series in the external momenta, instead of asymptotic 
expansions. In the Standard Model case, we confirm previous results derived on- 
shell and we obtain BR(B -> X s -f) = (3.62 ± 0.33) x 1CT 4 . In the case of the 
usual two-Higgs doublet model, we show that going from the leading to the next- 
to-leading order result improves the CLEO bound on the charged-Higgs mass from 
260 GeV to 380 GeV. This limit is very sensitive to the definition of errors and we 
carefully discuss the theoretical uncertainties. Finally, in the case of the two-Higgs 
doublet model in which both up- and down-type quarks couple to the same Higgs 
field, the theoretical prediction for BR(B — > X s ^) can be reduced by at most 20% 
with respect to the Standard Model value. 
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1 Introduction 



Interest in radiative B decays has been renewed in recent years, after the first mea- 
surements of the exclusive BR(B — > K*j) [|TJ and the inclusive BR(B —>■ X s 7) 
were performed by the CLEO collaboration. Theoretically these processes are inter- 
esting because they are sensitive to physics beyond the Fermi scale already at the 
leading order (LO) so that they are potentially well suited to probe new physics. 
Moreover these are "rare" processes with quite large branching ratios, of the or- 
der of 10 -4 , partly because they are strongly enhanced by the QCD corrections. 
Finally, the non-perturbative corrections to the inclusive decay turn out to be 
small || [|, D, |6], [7j, thus lending greater confidence on the perturbative calcula- 
tion. 

Experimentally, the measurement at the T(4S) @ 

BR(B -> X s7 ) = (2.32 ± 0.57 ± 0.35) x 10" 4 (1) 

was found in good agreement with the LO predictions || |J, which however suffer 
from a large theoretical uncertainty, ~ 30 — 40%, coming mainly from the choices 
of the renormalization and matching scales, which are not well defined at the LO. 

Recently the next-to-leading order (NLO) calculation in the Standard Model 
(SM) has been completed [1^, [11], |T^, |I3|| , raising the accuracy to the 10% level. 



The new predictions [0, Q turn out to be slightly larger than the CLEO mea- 
surement. However preliminary ALEPH results [15| now indicate a larger branching 
ratio measured at the Z° resonance, 

BR(B -> X s7 ) = (3.38 ± 0.74 ± 0.85) x 10" 4 . (2) 

Since this result is only preliminary and a complete analysis is now in progress, we 
will not use it in our discussion. However, if this result is confirmed, it shows that 
the experimental situation is not settled. 

In this paper, we present a new computation of the matching conditions of the 
AB = 1 magnetic and chromo-magnetic operators in the SM and in the two-Higgs- 
doublet model at the NLO, using an off-shell procedure. We confirm the previous 
results in the SM \TU, Tfi, [nj], obtained by on-shell matching, and present the charged 



Higgs contribution to the Wilson coefficients in the two-Higgs-doublet model. The 
paper is organised as follows. In sec. 0, we give some general formulae useful for im- 
plementing the matching conditions and discuss the peculiarities of our approach. 
Section |3| is devoted to some technical details of the calculation. In sec. |], the 
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relevant formulae for computing the NLO corrections to BR(B — > X s ^) are col- 
lected and a new phenomenological analysis in the context of the SM is presented. 
Sec. H contains the new coefficients computed in the two-Higgs doublet model and 
a discussion of the bounds enforced by the measurement of BR(B — > X s 7) on the 
parameter space (M#+,tan/3) of the model. Particular care is devoted to the def- 
inition of the theoretical errors, which play a major role in the phenomenological 
analysis, especially in the two-Higgs doublet model. Finally, in sec. ||, we summarise 
our results and conclusions. 



2 The Matching Conditions at the NLO 

In this section we discuss the determination of the matching condition of the AB = 1 
effective Hamiltonian H e ff at the NLO. 

The AB = 1 effective Hamiltonian is given by 

H eff = -^V*V tb Q(v)-C(ri , (3) 

where Gp is the Fermi coupling constant and V is the Cabibbo-Kobayashi-Maskawa 
(CKM) quark mixing matrix. Q(/i) is the row vector containing the renormalised 
dimension-six local operators which form the basis for the operator product expan- 
sion (OPE) of the AB = 1 weak-current product and C(p) are the corresponding 
Wilson coefficients. The determination of the effective Hamiltionan requires three 
steps: i) choice of the operator basis; ii) computation of the initial conditions for 
the Wilson coefficients by matching suitable amplitudes computed in the "full" and 
in the effective theory renormalised at a matching scale fi w ; iii) calculation of the 
relevant anomalous dimension matrix which determines the dependence of the Wil- 
son coefficients on the renormalization scale fi through the renormalization group 
equations. 

In the following we present the computation of the matching conditions for the 
(chromo-) magnetic operators in the case of the SM and of the two-Higgs-doublet 



model. The calculation in the SM, originally carried out in ref . [|HJ , has been already 



recently confirmed [16|, |14fl . The case of the two-Higgs doublet model is a novelty. 

A remarkable difference between our calculations and previous ones is that we 
work with the off-shell effective Hamiltonian, matching off-shell matrix elements, 



while the authors of refs. [fQ, 16, [L4] always use on-shell conditions. The advantage 
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of an off-shell matching relies on the fact that one can choose a suitable kinemati- 
cal configuration such that the relevant Feynman diagrams can be evaluated using 
ordinary Taylor expansions in the external momenta, in contrast to the on-shell 
case [16], in which asymptotic expansions [17] must be employed. However, this 



requires to work with an enlarged operator basis. 

The off-shell operator basis relevant for our calculation is [13, [1^, fLS 

Qi= {sLl^c L ){c L rt a b L ) 

Q2 = {sLlfiC L ){cLl^b L ) 

Qs= (s LlfM b L )E q (qrq) 
Qa = (s L J,t a b L ) E q (qYt a q) 

Qe= (s L ^ ulp t a b L )j: q (qrrYt a q) 

Q 7 = ^m^SLcr^bRjF^ 

(4) 

Qs = ^m b (s L a^t a b R )G^ K } 

Q 9 = j^s L pppb L 

Qio= ^s L {P,a^F^}b L 

Qn= -^s Llv b L D^ 

Q12 = ji^m b s L JP ipb R 

Qis= ^s L {p,a^t a G^}b L 

Qu= ^s L Yt a b L D»G% , 

where g s (e) is the strong (electromagnetic) coupling constant, qL,R are the chiral 
quark fields, G° v (F^) is the QCD (QED) field strength, is the SU(3) C x U(1) Q 
covariant derivative and t a are the colour matrices, normalised such that Tr (t a t b ^j = 
S ab /2. On-shell Q9,...,Qu are no longer independent operators and, by applying 
the equations of motion (EOM), they can be written as linear combinations of the 
operators Qi, Q 8 . 

In order to calculate the Wilson coefficients of the magnetic operators, one has 
to match a suitable amplitude computed both in the "full" and the effective the- 
ory between the same external states at a matching scale \x w . Schematically, the 
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condition is 

(^\J"4\b) = ( S1 \Q(n w )C(n w )\b) , (5) 

where J^Jl denotes the weak current product which is responsible for the AB = 1 
transitions in the "full" theory. 

At the NLO, retaining only the leading terms in 1/M^, the relevant amplitude 
computed in the "full" theory is given by 

(*y| J"Jl|6> = (si\Q {0) \b) ■ (r {0) + , (6) 

where (Q^) are the tree-level matrix elements of the operators in eq. (B). Intro- 
ducing the perturbative expansion 

C(ji) = (7(%) + ^(7«(/i) (T) 
and defining a matrix r such that 

(*r\$M\b) = (^ (0) + (si\Q {0) \b) , (8) 

the matching condition, eq. (H), becomes 

T (°) + {1) = r (0)T C (0) + ^i( r W T (7(°) + r (0)T (7 (1) ) (9) 
4tt 4tt 

namely 

(7(0) = f r (0)T\-lfQS) ^ = / (0)r|-l|j?(l) _ r (i)^( r (o)T)-iy(o)\ _ Qq\ 

Notice that i) in the presence of the magnetic operators, is not the identity 
matrix. Indeed the operator Q 2 has non-vanishing matrix element on Qu and Qn 
at the LO, through a one- loop diagram. Therefore the elements rfn and r^u 
do not vanish, ii) T"W contains terms proportional to log(m 2 /M^y) and similarly 
r^- 1 ' contains log(m 2 //i^), m 2 being the infrared regulator, so that the matching 
conditions depend on log(/i| / /M^ / ), but are independent of the infrared regulator. 
These log terms compensate a change of the matching scale [i w at the NLO, as 
stressed in ref. |14]]. Since the OPE must be valid on any matrix element, the 



Wilson coefficients cannot depend on the external states used to impose the matching 
conditions |2l|. For this reason our choice of matching on off-shell states is 
completely justified. 

According to eq. (|ToD, we first have to perform the tree- level matching by com- 
puting the relevant penguin diagrams for the various operator insertions. Only C 2 
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and C7-C14 are found to be non- vanishing at the LO. The LO matching of the off- 
shell effective Hamiltionian, (r^ T )~ 1 T ( -°\ has already been computed in ref. 
We have repeated the calculation and confirmed the results in eq. (27) of ref. 



Then we can proceed to calculate the finite parts of the two-loop diagrams, 
obtained by attaching gluons to the LO diagrams. Since we are only interested in 
on-shell magnetic operators, at the NLO we just need the coefficients of the operators 
which have a projection on the magnetic ones. Using the EOM, one finds 

C70S = C7 + ClQ 

^8 os = C*8 + Cl3 > (11) 

C 7os and C 8os being the coefficients of the on-shell magnetic and chromo- magnetic 
dipole operators respectively. Given the texture of r^°\ we have to compute only 
f/ 1} , and r§ } with i = 7, 8, 10, 13 and j = 1, . . . , 14. 

In the on-shell basis, it is more convenient to use the "effective" (chromo-) mag- 



netic coefficients, defined as || ||, [13 



C70S ~ C70S + ViCios 

4=1 
6 

^8 s = ^ 8 OS + ^2 Z iCi OS ' (12) 



4=1 



where y and z are scheme-dependent constant vectors, given by 

y Nnn = {QQ _l _1 _20 _80 

z» DR = (0,0,1,-^20,-^), (13) 
o 3 

in NDR. These definitions guarantee that the LO matrix element of the effective 
Hamiltonian for the transition b — > 57 (6 — > sg) is proportional to C%J (C*8os) on ly- 
Before outlining in the following section some technical details of the calculation, 
we discuss here some general features. 

1. We have regularized the ultraviolet divergences using naive dimensional regu- 
larization (NDR) with anti-commuting 75 and defined the renormalized oper- 
ators in the MS scheme. 



2. We use a background- field gauge [22| for QCD which avoids the appearance 



of non-gauge- invariant operators in the off-shell basis [23|. Moreover, in this 
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gauge, the external-gluon-field renormalization cancels out against the renor- 
malization of the strong-coupling constant which is present in the definition 
of the operators with an external gluon field. In the electroweak sector we 



employ the Fujikawa gauge |24|, B5| 



Flavour non-diagonal quark self-energies S^-(p) appear at one loop in the "full" 
theory. These can produce dimension-four operators in the effective theory, see 



e.g. ref. |TjJ. However it is well known [26] that the appearance of dimension- 



four operators can be avoided by choosing suitable renormalization conditions 
in the "full" theory, which cancel the momentum-independent part of the off- 
diagonal self-energies. This can be done by imposing that, for on-shell quark 
states q(p), gi(p)£y(p) = T, ij (p)q j (p) = 0, see e.g. ref. |f25|. 



4. Light-quark masses play the role of infrared regulators. However Feynman 
diagrams involving the three-gluon vertex present additional infrared diver- 
gences that are not regulated by this procedure. The simplest regulators in 
our approach are masses, thus we have computed these diagrams with a ficti- 
tious gluon mass. Since the gluon mass is bound to violate the QCD Ward- 
Takahashi identities at some order of the perturbative expansion [23|] , we have 



repeated the calculation regularising these infrared divergences using the di- 
mensional regularization. We have found that, in this particular calculation, 
the gluon mass is a "well-behaved" infrared regulator. 

5. A final observation concerns the matrix elements of Q3 and Q4. These op- 
erators have a scheme- dependent 0(oP s ) matrix element between h and 57, 
which is non- vanishing in NDR |27], B3]. Their contribution has to be explic- 



itly considered in the on-shell matching procedure [|1(], IB), |I4" |. This is however 
not necessarily the case for off-shell matching conditions, even if the same op- 
erators are also present in the off-shell basis. Indeed, these operators have 
vanishing off-shell initial conditions up to the NLO. The contribution coming 
from (^3,4 to the on-shell matching is provided off-shell by the 0(a s ) matrix 
element of Qii t n. The EOM account for the mismatch in the power of a s . 



3 Outline of the Calculation 

We work in an off-shell kinematical configuration in which all external momenta 
are assumed to be much smaller than the masses of the particles running inside the 



6 



loops. However, in this framework m s is kept only when it is needed as an infrared 
regulator. 

To extract C 7os (and analogously C 8os ) we notice that the most general vertex 
diagram having as external particles two fermions and one photon can be decom- 
posed in 12 x 2 objects, where 12 are the possible combinations of one Lorentz index 
and two momenta (like, for example, p with p and q independent momenta) while 
2 takes into account the two possible chiralities. Using standard techniques |29| 
one can construct projectors for each of these 24 objects and then take the appro- 
priate combination of them which, after the use of EOM and the Gordon identity, 
contributes to the on-shell magnetic operator. We observe that these projectors 
depend on the kinematical configuration assumed. We find it convenient to work 
in a "scaled" configuration defined by q 2 = 0, p ■ q = p 2 /2 where p and q are the 
momenta carried by the b quark and the photon respectively, under the assumption 
p 2 <C m 2 , ml. In this configuration the projector for the tensor part of the on-shell 
(chromo-) magnetic operator is 

^7,8) = 2 p4 ft _ 2) {P 2 -ft -P 2 l^-m b t*# + 2ft jf -2p 2 p» 

+ rribnpq^ — m^n^q^ + n0q^j (14) 

and works by contracting it with the amplitudes and taking the trace. In eq. ( |14"D 
n is the dimension of the space-time and p is the index carried by the photon. As 
it can be seen from eq. fll4l) the projector contains powers of p 2 in the denominator. 
However, because of our assumption p 2 <C m 2 , m 2 , we can perform ordinary Taylor 
series to eliminate the unphysical poles when p 2 — » 0. Consequently, the actual 
computation of the graphs both for the "full" and the effective theory is reduced to 
the calculation of massive self-energy integrals with no external momenta. Explicit 
results on two-loop integrals of this kind are known [pOf| . 

The use of an off-shell configuration induces in the result for the "full" theory 
the appearance of terms that behave as l/m 2 , 1/m 2 , log m 2 , log as m s , nij ^ 0. 
However, the effective theory has the same infrared behaviour and there is a one- 
to-one correspondence between m s , — > divergent contributions in the diagrams 
of the "full" and those of the effective theory. It is worth reminding that this cor- 
respondence works in the off-shell operator basis of eq. @. In particular, diagrams 
like (la) in fig. [I] have a counterpart in the effective theory, depicted in (lb), that 
requires the introduction of a b — > sg^ vertex. Such a vertex is indeed present in 
the operators Qg, Qio, Q12, Q13 while it does not appear in the on-shell operator 
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Figure 1: Examples of the correspondence between diagrams in the "full" theory and in 
the off-shell effective theory. 

basis. The corresponding contribution can be worked out either through the use 
of the Ward identities that relate the 1PI b — ► sj, b —>■ sg vertex to the b — > sg^ 
one |3l] , or by direct computation of the Feynman rules from the relevant operators. 
In an analogous way, topology (2a) requires a b — > s transition vertex in the effective 
theory. The leading behaviour of (2a) as m s ,mb — > is 0(l/m%), 0(l/mf) and this 
is eliminated if the "full" theory is renormalized on-shell ]25| (see discussion in point 
3 of section 2). Alternatively, if the vanishing of the off-diagonal mass terms is not 
enforced, then the introduction of the dimension-four operators 

M 1= jt*s L pb L , 

, _ (15) 
M 2 = m b s L b R , 

in the diagram (2b) will provide such a cancellation. However, topology (2a) also 
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shows an infrared logarithmic behaviour in the light quark masses. This cannot be 
eliminated through renormalization but it is actually cancelled when the operators 
Qg and Qu are inserted in diagram (2b). 

It should be noticed that the cancellation of infrared m s , rrib contributions be- 
tween "full" and effective theory works before enforcing the GIM mechanism. There- 
fore the charm diagrams do not present any additional difficulty The charm con- 
tribution to C*7 og can be obtained by replacing m t with m c and taking the limit 
m c — > 0. After performing this operation one is left with some log m c terms. How- 
ever, the effective theory contains also the contribution of the operator Q2 that 
provides the cancellation of these logarithmic terms through two-loop diagrams. 

As already said, the procedure for extracting C$ os is completely analogous, the 
only difference being that instead of introducing a b —>■ sg^f vertex, one needs a 
b — > sgg vertex. This is present in the operators Qg, Q12, Q13, Qu- Again, 
the infrared divergent contributions associated with the gluon cancel diagram by 
diagram between the "full" and the effective theory. 



4 Standard Model Results 



In this section we present the result for the SM. The analysis is performed in the on- 
shell basis obtained by applying the EOM to the operators Qg— Qu- All the Wilson 
coefficients are now intended in the on-shell basis, but for simplicity we drop the 
subscript OS. The NLO expression for the branching ratio of the inclusive radiative 
decay B 



A s 7 has been presented in ref. [[13 
BR^ = BR(B -> A s7 ) = BR(B -»■ 

x (\D\ 2 + A) I 1 



X r eu e 



v: s v tb 2 



V c 

sNP 
°SL 



cb 



irf(z)K(z) m 2 b 



+ 



8. 



NP 



rnr, 



+ 



s 



(16) 

% m b m c J 

Here and throughout the paper m q (q = c, b, t) denotes the quark pole mass, while 



m„ denotes the MS running mass, 



m q {m q ) 



_ A a a (m q ) 

3 7T 



(17) 



Also z = m 2 c jm\ and f(z) and k(z) are the phase-space factor and the QCD correc- 

2a s (fi b ) h(z) 



tion |32| for the semileptonic decay, 

k(z) = 1 



37T f( Z ) 



18) 
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f(z) = l-8z + 8z 3 - z 4 -12z 2 \nz 



(19) 



h(z) 



., ,'25 239 25 2 \ , / 4 2 17 ^ 

~ ^) ( 4" ~ IT* + I J + [ + ~ 3 + ~3~ , 



+z 2 \n 2 z (36 + z 2 ) + (1 - 



17 64 17 o , , , 



-4(1 + 30^ 2 + ^ 4 ) In z ln(l - z) - (1 + 16z 2 + z 4 ) 6Li 2 (z) - tt 

1 - \ 



-32z 3/2 (l + 



7T 



4Li 2 (v^) +4Li 2 (-v/?) - 21nz In 




(20) 



In eq. ( pf) /4 is the renormalization scale relevant to the semileptonic process, a 
priori different than the renormalization scale relevant to the radiative decay, as 
stressed in ref. 

In order to relate the quark decay rate to the actual hadronic process, we have 
included the corrections obtained with the method of the Heavy-Quark Effective 



Theory (HQET), for a review, see e.g. ref. p3 |. In terms of the contributions to a 
heavy hadron mass from kinetic energy (Ai) and chromo-magnetic interactions (A 2 ), 
the 1/m 2 corrections to the semileptonic and radiative decays are 



?NP Ai 3 

d SL = TT + o A 2 



[i-zY 



Ai 9. 

— — — Ao 
2 2 



(21) 



(22) 



Notice that the dependence on the poorly-known parameter Ai drops out in eq. 



We take A 5 



m 



B" 



-m 2 B )/A = 0.12 GeV 2 , as extracted from B-meson mass splitting. 



Voloshin H has considered also 1/m 2 non-perturbative corrections coming from 
interference between the operators 2 and Of. With the correct sign found in 
ref. M, these are given by 



rJVP 



1 A 2 



(o) 



d 0) 



(o)' 



(23) 



where C 2 °' ) and Cf^ are the corresponding Wilson coefficients at the LO computed 
at the scale It appears that higher-order terms in 1/m 2 are suppressed 0, and 
therefore non-perturbative corrections are presumably under control. 

In eq. (^), A is the correction coming from the bremsstrahlung process b — > s'yg, 



with a cutoff on the photon energy E~ > (1 — 5)rrib/2 [|TT], ^4], 
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.4 



^ e -a s (Mb)ln<5(7+21n5)/37r _ -A | Q (0)eff ^ 1 2 
+ «M £ C (0 e// ( ^ )c (0) .// (A<6)/ ((5) . 



(24) 



The coefficients fij(5), which have been computed in ref. JTT], 34], are defined as in 
ref. ||13|| . To conform to previous literature, we define a "total" inclusive decay by 
choosing 5 = 0.99, although for comparison we will also present results for 5 = 0.90. 
The numerical values of the coefficients fij(5), for 5 = 0.99 and y/z = 0.29, areQ 



0.0009 



fij — 



-0.0113 
0.0340 



-0.0035 
0.0210 
0.0140 



0.0006 


-0.0459 


-0.0600 


-0.0035 


0.2754 


0.3599 


-0.0047 


0.3277 


0.0666 


0.0088 


-0.0546 


0.1570 




1.9369 


0.9506 






0.9305 



-0.0030 
0.0182 
0.0421 

-0.0070 
0.5926 
0.0002 
3.4211 



0.0010 \ 
-0.0061 
-0.0140 

0.0023 
-0.1975 
-0.0001 

0.3897 

3.2162 J 



(25) 



Finally 



D = Cf e// (/i b ) + 



An 



(0)e//, 



(26) 



The coefficients r», obtained from an NLO matrix element calculation, have been 
calculated in ref. |l!|] for i — 1, 2, 7, 8. They can be found in ref. fl3| , together with 



the corresponding values of "f$ e ^ ' ■ The dependence on the renormalization 
scale relevant for the B —>■ Xg'j process, cancels out in D at the leading order in 
a s . The relations between the Wilson coefficients at the low-energy scale \x\, and the 
high-energy scale fi w , obtained by renormalization group evolution, can be found in 
refs. []13| , |14| . We are interested here in the expressions of the Wilson coefficients at 
the scale fi w , at which the "full" theory is matched into an effective theory with five 
quark flavours. In the SM model, they are given by 



(0)e//, 



4tt 



-C 



(27) 



1 We thank M. Misiak and N. Pott for pointing out to us misprints in the expressions of /12 in 
ref. Ill and / 78 in ref. M. 



11 



(0), 




1 



for i = 1,3, 4, 5, 6 
for i — 2 
for i = 7,8, 



F 7 (1) (x) 



F 8 (1) (x) 



a; (7 - 5a; - 8a; 2 ) x 2 (3x - 2) , 

+ — 77" 7T77- In x 



24(x-l) 3 4(x-l) 4 
x(2 + 5x — x 2 ) 



3x 2 



m 2 (fi w ) 



In x 



M 2 



(28) 

(29) 
(30) 
(31) 



The NLO top-quark running mass at the scale fi w is given by 



m t {Hw) = m t (m t ) 





~<0 


a s (fi w ) 


2/3 


_a s (m t )_ 





4tt 2^o \7o m A)J V«.N 



A, = 11 --71,, ft = 102-^71, 



(32) 
(33) 

(34) 

3 9 J ' V ; 

and it is related to the pole mass by eq. fllTD - Here n, is the effective number of 
quark flavours. We have computed the NLO corrections using the method illustrated 



7o m 



7i m 



404 40 



-n, 



in sect. 3 and confirmed the results of ref. flit) 



15 + 6 In 







E(x) 



for i = 1 
for i = 2, 3, 5, 6 
for z = 4 



2 , 2 j n % 

3 ^ 3 

Gi(a;) + A<(x)lii^f- for z = 7 



(35) 



- . x -18 + 11a; + a; 2 a; 2 15 - 16a; + 4x 2 , 2, 

-E(x) = — -7 H 1 : 77 In a; In a;. 

V ; 12(a;-l) 3 6(x-l) 4 3 



(36) 



G 7 (x) 



■16a; 4 - 122x 3 + 80x 2 - 8 



x T / 1\ 6x 4 + 46a; 3 - 28a; 2 , 2 

L n 1 —) + 3 (x-i)» lni 



+ 



9(x- l) 4 

-102a; 5 - 588a; 4 - 2262a; 3 + 3244a; 2 - 1364a; + 208 



}l(a; - l) 5 



In a; 



2 We are grateful to F. Borzumati and C. Greub for pointing out to us the omission of the 
logarithmic corrections to C^\ e ^ {Hw) m the preprint version of this work (see also ref. p5|). 
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Mb 
[GeV] 


[GeV] 


[GeV] 


m t 
[GeV] 




m c /m b 


m b - m c 
[GeV] 


BRsl 


a s {M z ) 
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4.8 
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Central values 
0.95 I 0.29 


3.39 


0.1049 


0.118 


130.3 


+4.8 

-2.4 


+4.8 
-2.4 


+80.33 
-40.16 


±5 


Errors/Ranges 
±0.03 | ±0.02 


±0.04 


±0.0046 


±0.003 


±2.3 


-190/ 
+24 /o 




+9.3o/ 
-9.7 /o 


±1.6% 


Uncertainty in BR~ f 
±3.2% | tf 5 3 % 


(LO) 

-0.07q/ 
+0.08 70 


±4.4% 


+2.30/ 
-2.2 /o 


-1.70/ 
+ 1.8 /o 


+°-°% 
-3.5 /0 


-1.4 
+2.1 


-Uw 
+2.1 /o 


±1.4% 


Jncertaint 
±3.2% 


y in BR 1 

+4.50/ 

-3.8 /o 


(NLO) 
TO.4% 


±4.4% 


+2.90/ 
-2.8 /o 


+1.8% 



Table 1: The input parameters entering the calculation of BRj, their central values, 
their errors, their influence on the uncertainty in BR^ at the LO and NLO. The other 
constants entering the calculation are M w = 80.33 GeV, A2 = 0.12 GeV 2 , m s /mb = 50, 
and S = 0.99. 



1646a; 4 + 12205a; 3 - 10740a; 2 + 2509a; - 436 

H 7 vi (37) 

486(a; - l) 4 K ' 

„ . . -4a; 4 + 40a; 3 + 41a; 2 + x T ( 1 \ -17a; 3 - 31a; 2 , „ 

G8W = S^iji L H 1 -xJ + 2 («-iy * 

-210a; 5 + 1086a; 4 + 4893a; 3 + 2857a; 2 - 1994a; + 280 , 

H ; r~z mX 

216(3;- l) 5 

737a; 4 - 14102a; 3 - 28209x 2 + 610a; - 508 
+ 1296(x-l) 4 (38) 

a / \ 208 - 1111a; + 1086a; 2 + 383a; 3 + 82x 4 2x 2 (14 - 23x - 3x 2 ) , 

A 7 x = — — 4 + K — — 5 J -\nx 39 

81(x — l) 4 3(x — 1) & 

* / \ 140 - 902a; - 1509a; 2 - 398a; 3 + 77a; 4 a; 2 (31 + 17 x) , 

A$(x) = -, r-j 1 rrr^lnx (40) 

V 1 108(x-l) 4 2(x-l) 5 V ' 

The ji w dependence in all terms of order a s in eq. (^) cancels out because eqs. ( |3~9"D 
and (HOT) satisfy the relation 



A - = ^^t + \p*'" c ?' («) 



The values of the parameters chosen in our numerical analysis are shown in 
table 1. The value of a s (Mz), obtained from an average of low- and high-energy 
measurements, is taken from ref. [BH]. From this we calculate a s at a generic scale 
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\x with the NLO result 

aM = «jmU-^JMA^) „ = 1+A) ^ ln JL (42 ) 

v \ p 47r v J 2it M z 

The top quark mass is obtained by combining the results by CDF [|37j (m t = 
176.8 ± 4.4 stat ± 4.8 syst GeV) and DO g| (m t = 173.3 ± 5.6 sia * ± 6.2 syst GeV). The 
difference between the bottom and charm mass is obtained from the HQET relation 
between the quark and hadron masses, m b — m c = 3.39±0.03±0.03 GeV []55[ ], where 
the first error is from the uncertainty in Ai and the second one from higher-order 
corrections. The ratio m c /mb is chosen in the range 0.29 ± 0.02, which corresponds 
to m b = 4.77 ±0.15 GeV. 

The uncertainty on the electromagnetic coupling a e corresponds to the arbitrari- 
ness of choosing the renormalization scale between rrib and M w , as QED corrections 
have not been included. 

The experimental situation on the semileptonic branching ratio BRsl is rather 
controversial (for a review, see ref. [(KJ). The two measurements at the T(4S) give 
BR SL = (10.49 ± 0.17 ± 0.43)% (CLEO g|) and BR SL = (9.7 ± 0.5 ± 0.4)% (AR- 
GUS |52|]). On the other hand, the average of the four LEP experiments gives f4*3] , fill 
BR Z SL = (10.95 ± 0.13 ± 0.29)%. However the B samples at the Z° also contain 
B hadrons which have typically shorter lifetimes than B mesons. Following Neu- 
bert P9f, we can correct BR§ L using the average lifetime of the corresponding B 



mixture and obtain BRsl = (11.23±0.34)% at Z°. The low- and high-energy results 
show a certain discrepancy. However, a reanalysis of the LEP results is finding new 



systematic errors and a new dedicated analysis is in progress ||44|| . For this reason 
we prefer to use only the CLEO result. 

In terms of the Wolfenstein parameters A, p, and r] related to the standard 
parametrisation [[36] by 

si2 = A, = \(p-i v ) (43) 

S23 

and using unitarity, the relevant CKM factor can be written as 



V t * s V th 2 



cb 



1 + (2p - 1)A 2 . (44) 



For A = 0.221 and p in the range between —0.3 and 0.3, we obtain the value shown in 
table 1. The new physics contributions that we will consider in sect. 5 can modify the 
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values of p obtained from fits of different flavour-physics measurements. However, 
since we allow a rather large range for p, these effects are included in the error. 

To estimate higher-order effects, we vary the different renormalization scales 
Pb, and p w in the ranges given in table 1, which correspond to shifts in the logarithms 
of the scale by an amount ± In 2 around their central values Pb = Pb = m-b, AV = M w . 
In the NLO calculation we have dropped all explicit 0{a 2 s ) terms in eq. (|16|). To allow 
a better comparison, we also perform the calculation in LO approximation, where 
0(a s ) terms are set to zero and a s is evolved using the one-loop renormalization 
group equation. 

The result of the numerical evaluation of eq. (|16|) is 

BRy = (2.66 ± 0.74) x 1CT 4 (LO) (45) 



BRy = (3.62 ± 0.33) x 10" 4 (NLO, 5 = 0.99) 

(3.54 ± 0.32) x 10~ 4 (NLO, 5 = 0.90) (46) 

The error is obtained by adding in quadratures the errors due to the different input 
parameters shown in table 1. If the error is asymmetric, we always choose the 
larger uncertainty. The pb scale dependence dominates the error at the LO, but 
it is drastically reduced at the NLO. This is illustrated in fig. 2. The result in 



eq. ( [46]) is in agreement with those presented in refs. ]1^, [L4J]. The slightly different 
Pb dependence we are finding comes from our choice of defining fhb in eq. ( P"6"[) at 
the scale pb (where the operator Q 7 is defined), instead of m b as in refs. |L3|, [14 



Of course this just corresponds to a higher-order effect, and this ambiguity is taken 
into account by the theoretical uncertainty. 

Given the controversial situation in BRsl, it is useful to quote also the prediction 
for the ratio BR 7 /BR SL (for 5 = 0.99) 

BRJBRsl = (2.54 ± 0.70) x 10" 3 (LO) (47) 

BR 1 /BR SL = (3.45 ± 0.27) x 10^ 3 (NLO) (48) 



5 Two-Higgs Doublet Models 

In this section we investigate the prediction for BR y in models with two Higgs 
doublets. The new ingredient with respect to the SM is the presence of a charged 
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Figure 2: The fi b scale dependence of BR y (5 = 0.99) at LO (lower curves) and NLO 
(upper curves). The solid line represents the case in which all input parameters are 
fixed to their central values. The dashed lines represent the errors coming from all input 
parameters, other than fx b , added in quadratures. 

Higgs scalar particle which gives new contributions to the Wilson coefficients of 
the effective theory. The interaction between quarks and the charged Higgs if ± is 
defined by the Lagrangian 

£ = (2V2G F ) 1/2 J2 «< Uum Ui Vi,^-^ - A d V l3 m di ^-^-\ d 3 H + + h.c. (49) 

i,j=l \ - 

Here i,j are generation indices, m u ^ are quark masses, and V is the CKM matrix. 
We are considering only models in which FCNCs vanish at tree level ||5] . Therefore 
we have two possibilities which are parametrised by the mo del- dependent coefficients 
A u and A d . The first possibility, referred to as Model I, is that both up and down 
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quarks get their masses from the same Higgs doublet H 2 , and 

A u = A d = l/tan/3 , (50) 

where tan (3 = (H®) / (H®) . The case of Model II is that up quarks get their masses 
from Yukawa couplings to H 2 , while down quarks get masses from couplings to H 1 , 
and 

A u = -1/A d = 1/ tan/3 . (51) 

The latter case is particularly interesting because it is realized in the minimal su- 
persymmetric extension of the SM. 

In the language of the effective theory, it is easy to include the charged-Higgs 
effect. We just need to compute the new Wilson coefficients at the matching scale 
/i w , since they contain all the relevant ultraviolet information. At the LO, the 
charged-Higgs contribution, to be added to eq. ( |28| ) is given by f46| 



SCl 0)eff (ji w ) = Q z = l,...,6 (52) 

SCit ff (fi w ) = - A u A d F%{y) (53) 

v 2/(3 - 5y) y(3y - 2) 



12(2/ - l) 2 6(y -r 



^) = f^-^ln, (55) 

»=*^, (56) 

and F$ are defined in eqs. (H)-©. 

At the NLO, the charged-Higgs contributions to the Wilson coefficients are 

<JCf Mf (nw) = i = 1,2,3,5,6 (57) 
5Cf )e// (M =£"(?/) (58) 

^ 1)£// (^) = G?(y) + Af( y )ln||- V(y) 

^ 1)£// (^) = Gf(y) + Af(y)lB^-^(y) (59) 
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with 



A A 4 

A d A u -y 



4(-3 + 7y-2y2) / i\ 8-Uy-3y\ 2 



2(-3-y + 12y 2 -2y 3 ) 1 7-13y + 2y 2 
H — — ln y + 



3(y-l) 4 

y(18-37y + 8y 2 ) T . 



(y - 1) ; 



(y-i) z 



Li 2 1 



y, 



+ 



y(-14 + 23y + 3y 2 ) 

(y - 1) 5 



-50 + 251y - 174y 2 - 192y 3 + 21y 4 , 

+ w^f h,v 

797 - 5436y + 7569y 2 - 1202y 3 
+ 108(j/ - l) 4 , 



In 2 y 



(60) 



Af(y) 



A d A u -y 



+Al-y 



21-47y + 8y 2 , 2(-8 + 1% + 3y 2 ) 



ln y 



(y - i) 3 (y - 1) 4 

-31 - 18y + 135y 2 - 1 V y(U - 23y - 3y 2 ) 



6(y-l) 4 



+ 



(y-if 



lny 



(61) 



G?(v) 



A d A u -y 



-36 + 2%-lV) / 1\ 19 + 171/ , 2 



+ -3- 187y + 12y 2 - 14y 3 ) ^ + 3(143 - 4% + 29y 2 ) 



4(2/- l) 4 
y(30-17y + 13y 2 ) 



8(1/ -l) 3 



+Al-y 



(y-i) z 



Li, 1 - 



y, 



2/(31 + I7y) 2 



-226 + 817y + 1353y 2 + 318y 3 + 42y 4 , 

+ m^Ty' ln9 

1130 - 181531/ + 7650i/ 2 - 4451y 3 



216(i/ -1)^ 



(62) 



Af(y) 



A d A u -y 



11 - 16y + 7y 2 _ 19 + 17y 
2(1/ -l) 3 (y-i) 4 



-38 - 261y + 18y 2 - 7y 3 + y(31 + 17y) ^ 



6(y - l) 4 



(y-i) E 



2/(16 - 29y + 7y 2 ) + y(3y - 2) ^ ^ 



36(y - 1) J 



6(y - l) 4 



(63) 
(64) 
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M H [GeV] 

Figure 3: Predicted BR^(5 = 0.99) in Model II as a function of the charged-Higgs 
mass Mjj, for tan/3 = 2 at LO (solid line). The dashed lines represent the uncertainty 
obtained by adding in quadratures the errors coming from the different input parameters. 
The dot-dashed lines define the uncertainty range with the "scanning" method (see text). 

The fi w dependence cancels in terms of 0(a> s ), since the contribution satisfies 
eq. ©. 

In figs. [|and [|we show the predicted value of BR^ in Model II as a function of the 
charged-Higgs mass Mh, for tan/3 = 2 at LO and NLO with 5 = 0.99. As apparent 
from these figures, the decoupling at large Mh occurs very slowly. Moreover, the 
contribution always increases the SM value of BR^. Since the SM prediction 
is already quite higher than the CLEO measurement, one can obtain very stringent 
bounds on Mh ||47|| . These bounds are very sensitive on the errors of the theoretical 
prediction. Improving the calculation to the NLO has important effects on these 
bounds, since the theoretical error is significantly reduced, as shown in figs. 3 and 
4. 
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M H [GeV] 

Figure 4: Predicted BR^(5 = 0.99) in Model II as a function of the charged-Higgs 
mass M H , for tan/5 = 2 at NLO (solid line). The dashed lines represent the uncertainty 
obtained by adding in quadratures the errors coming from the different input parameters. 
The dot-dashed lines define the uncertainty range with the "scanning" method (see text). 

Because of the large sensitivity on errors, we will use two different methods to 
define the bounds on M H . We define a "Gaussian" standard deviation by adding 
different input errors in quadratures. We also define a "scanning" method of estimat- 
ing errors by varying independently all input parameters within the ranges defined 
in table 1. Figures 3 and 4 show the errors on the theoretical prediction using 
one "Gaussian" standard deviation (dashed lines), or using the "scanning" method 
(dot-dashed lines). If we apply the "scanning" method to the SM, we obtain (for 
5 = 0.99) 

BR 1 = (2.66ll$j) x !0~ 4 (LO) (65) 
BR, = (3.62lg$) x 10~ 4 (NLO) (66) 
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M H [GeV] 

Figure 5: Exclusion regions in the plane BR^j BR SL versus M H for tan/3 = 1 (dashed 
lines), 2 (solid lines), 5 (dot-dashed lines) obtained by our NLO calculation with 5 = 
0.90. The upper lines correspond to a theoretical uncertainty defined by one "Gaussian" 
standard deviation, and the lower lines to the "scanning" method (see text). 

BRJBRsl = (2.54+^) x lO" 3 (LO) (67) 

BRy/BRsL = (3.45+°;^) x 10~ 3 (NLO) (68) 

This has to be compared with the "Gaussian" standard deviation error given in 
eqs. (HP-fl). 

Next, we exclude a value of M# if the minimum of BR^ inside the "Gaussian" 
(or "scanning") band lies above the present CLEO 95% CL upper bound BR^ < 
4.2 x 10~ 4 0. The corresponding limits on M# are shown in table 2, for different 
values of tan (3. If we restrict our consideration to the case tan (3 > 1, the dependence 
on tan/3 is very mild, and the limit practically saturates for tan/3 ^ 2. Notice 
that, for a "Gaussian" error, going from the LO to the NLO calculation allows 
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tan p 


theor. error 


M H (LO) 


71 JT /~\TT / ~\ \ 

Mh (NLO) 


1 /T /ATT / A \ 

M H (NLO) 






[GeVJ 


[■/-"i t rl / r r\ r\r\\ 

[GeVJ (5 = 0.99) 


r s~~i i rl / C" n aa\ 

[GeVJ (6 = 0.90) 


1 


Gaussian 




A 70 

473 


a r\ o 

408 


2 


"Gaussian" 


268 


440 


377 


5 


"Gaussian" 


258 


430 


368 


1 


"scanning" 


240 


316 


276 


2 


"scanning" 


208 


287 


249 


5 


"scanning" 


198 


278 


242 



Table 2: Lower bounds on the charged-Higgs mass Mh from the CLEO 95% CL limit 
BR 7 < 4.2 x 10 -4 , using LO or NLO calculations. The theoretical uncertainty is 
estimated either with one "Gaussian" standard deviation, with input parameter errors 
added in quadratures, or with the "scanning" method, by varying independently the 
input parameters within their range. The NLO results are presented for two different 
values of the cut on the photon energy 5. 



an improvement of the limit on Mh from about 260 GeV to about 370 GeV, for 
large tan/5, even in the more conservative case of 5 = 0.90. Had we neglected 
the NLO Higgs matching condition, but kept the remaining NLO corrections (as 
done in ref. |f48|), we would find a much stronger limit M# > 460 GeV (tan/3 = 5, 
"Gaussian" error). 

It may seem disturbing to integrate out the top quark, the W, and a charged 
Higgs as heavy as 300 GeV or more at the same scale fi w . However, since we have 
completed the NLO calculation, the residual dependence on \i w is very small. A 
resummation of the leading logarithms of Mh/ttiw has been performed in ref. |49|. 

Certainly one of the reasons we are able to obtain such strong bounds on Mh 
is the poor agreement between SM prediction and CLEO measurement. In view of 
the preliminary ALEPH result on BR^ and the controversial situation on BRsl, 
it is more appropriate to present our results in a fashion which is independent of 
the present measurements. This is done in fig. 5, which shows the bound on Mh 
according to the "Gaussian" and "scanning" methods, for any given value of the 
experimental upper bound on BR^j BRsl in the more conservative case 5 = 0.90. 
The region below the curves is excluded by our theoretical calculation. As future 
measurements on BR^j BRsl improve, on can simply read off the limit on Mh from 
the lines of fig. 5. 

It should be remarked that, although Model II corresponds to the supersymmet- 
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Figure 6: The ratio between BR y (5 = 0.90) predicted in Model I and the SM value, for 
the minimum tan/3 allowed by the Rb measurement (solid line), for tan/3 = 2 (dashed 
line) and tan/3 = 3 (dot-dashed line). The solid line is interrupted at tan/3 = 1. 



ric Higgs sector, these bounds do not directly apply to the supersymmetric extension 
of the SM. The reason is that contributions coming from other new particles can 
partly compensate the effect of the charged Higgs. In particular, in the supersym- 
metric limit, there is an exact cancellation of different contributions 1501. In the 



realistic case of broken supersymmetry, this cancellation is spoiled but, if charginos 
and stops are light, it may still be partially effective. Nevertheless, the bound pre- 
sented here severely restricts the allowed region of parameter space. 

In the case of Model I, the charged-Higgs contribution reduces the value of BR^ 
and therefore we cannot obtain any significant bound. On the other hand, it is 
interesting that new physics effects can bring the prediction for BR 1 closer to the 
CLEO value. 

A significant effect can only be expected for small tan/3, since in Model I all 
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charged- Higgs contributions vanish in the limit of large tan /3, as tan -2 (3. When 
tan (3 is reduced, the top Yukawa coupling grows, and one eventually is limited by 
new contributions to B°-B° mixing, and occurrence of a Landau pole at interme- 
diate energies (see e.g. ref. [51]). We have checked that the strongest experimental 
lower bound on tan/3 comes from R b , the ratio between the bottom and hadronic 
Z° widths^. Using the 95% CL lower limit > 0.2155 obtained from a combined 
LEP+SLD analysis (R b = 0.2170 ± 0.009 @) and the smallest top-quark mass 
value inside our range (m t = 170 GeV), we derive a lower bound on tan/3 of 1.8, 
1.4, 1.0 for Mh = 85, 200, and 425 GeV, respectively. For these values of tan/? 
and M H , the B°-B° mixing parameter is increased by about 30% with respect to 
the SM. Such an increase is compatible with experiments, given that the theoretical 
uncertainty on the hadronic parameter F^Bb is at least 40%. In correspondence 
to this minimum value of tan/? we obtain the maximum reduction factor of BR^ 
with respect to the SM value, plotted in fig. 6. This reduction can be at best about 
20%. For comparison we also show the theoretical prediction on BR 7 / (BR^sm for 
different values of tan (3. The lines in fig. 6 correspond to the central values of the 
input parameters defined in table 1, and error effects have not been included. 



6 Conclusions 

We have presented the 0(a s ) QCD corrections to the matching conditions of the 
flavour-violating magnetic and chromo-magnetic operators both in the SM and in 
two-Higgs doublet model. We have worked with an off-shell effective Hamiltonian, 
and this has allowed us to compute the relevant Feynman diagrams without the use 
of asymptotic expansions. 

Following this different method we have confirmed the results for the matching 



conditions in the SM first found in ref. |fL0|1 , and recently recomputed in refs. [JIB], [14 
Our numerical analysis for the branching ratio of B — > X s 7 gives 

BRy = (3.62 ± 0.33) x 10" 3 (69) 

in agreement with previous literature [|T^, |H]j . 

We have then extended our results to the case of two-Higgs doublet models. In 
the case of Model II, the charged-Higgs contribution always increases the SM pre- 
diction. Since the present CLEO measurement lies below the SM result, B — > X s ^y 



3 We thank C. Wagner for providing us with a numerical calculation of Rb- 
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provides quite stringent bounds on the charged-Higgs mass Mu- These bounds are 
very sensitive to the theoretical errors. The improvement of the calculation to the 
NLO significantly reduces the theoretical uncertainty and therefore tightens con- 
siderably the Mu bounds. Our results are summarised in fig. 5, which shows the 
region in the plane BR^j BR SL versus M H excluded by our theoretical calculation. 
An experimental upper bound on BR^/BRsl can then be easily translated into a 
bound on Mu- For instance, using the present CLEO experimental results and a 
theoretical uncertainty defined by one "Gaussian" standard deviation, we find that 
charged-Higgs masses below 370 GeV are excluded. This result improves the bound 
of about 260 GeV obtained by LO calculations. We recall that the experimental 
results for BR$l show a certain discrepancy between low- and high-energy mea- 
surements (as discussed in sect. 4) and the preliminary ALEPH result for B — > 
is not in perfect agreement with the CLEO measurements. As shown in fig. 5, de- 
velopments in the experimental situation can have an important effect on the Mu 
bounds. 

This bound on Mu does not directly apply to the case of supersymmetric models, 
since contributions from other new particles can compensate the charged Higgs 
effect, and this is indeed known to happen in certain cases [[50]. Nevertheless, it is 



clear that this bound is so strong that it severely restricts the allowed parameter 
space. This has important consequences on the prediction of new particle masses 
and on their observability at future colliders. However, it has to be remarked once 
again that the impressive restrictiveness of our bound on Mu is partly a consequence 
of the poor agreement between theory and the CLEO measurement. Developments 
on the experimental side or new understanding of non-perturbative effects may then 
partially relax this bound. 

Finally we have investigated the Model I version of theories with two-Higgs 
doublets. In this case, the charged-Higgs contribution reduces the SM value for 
BR^. Therefore one cannot derive useful bounds on Mu , but instead the theoretical 
prediction can be brought closer to the CLEO measurement. However, in view of 
the limits on R b , the Model I charged Higgs can reduce the SM prediction for BR^ 
by at most 20%. 

We wish to thank F. Borzumati, P. Chankowski, K. Chetyrkin, C. Greub, 
M. Mangano, G. Martinelli, M. Misiak, M. Neubert, M. Testa, and C. Wagner 
for useful conversations and D. Schlatter and H. Wachsmuth for discussions of the 
experimental results. We also acknowledge useful communications with A. Buras, 
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G. Isidori, and N. Pott. One of us (G.D.) would like to thank the Physics Department 
of New York University for the kind ospitality during the initial stage of this project. 



Note added: The calculation of the matching condition of Qj for the two-Higgs 



doublet model II, in the large tan (3 limit, has also been recently presented in ref. |52 



and the result agrees with our eq. (|60|). In the final version of this paper, various 
formulae entering the calculation of BR(B — > X s ^) have been corrected, following 
suggestions by F. Borzumati, C. Greub, M. Misiak and M. Neubert. 
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